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ABSTRACT 

Aims. We calculate Keplerian (mass shedding) configurations of rigidly rotating neutron stars and quark stars with 
crusts. We check the validity of empirical formula for Keplerian frequency, /k, proposed by Lattimer & Prakash, 
f K {M) = C (M/M Q ) 1/2 (R/10 km) _3/2 , where M is the (gravitational) mass of Keplerian configuration, R is the 
(circumferential) radius of the non-rotating configuration of the same gravitational mass, and C = 1.04 kHz. 
Methods. Numerical calculations are performed using precise 2-D codes based on the multi-domain spectral methods. 
We use a representative set of equations of state (EOSs) of neutron stars and quark stars. 

Results. We show that the empirical formula for fn{M) holds within a few percent for neutron stars with realistic 
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for an EOS, and C = Cns = 1-08 kHz. Similar precision is obtained for quark stars with 0.5 Mq < M < 0.9 M„ 
For maximal crust masses we obtain Cqs = 1.15 kHz, and the value of Cqs is not very sensitive to the crust mass. 
All our C's are significantly larger than the analytic value from the relativistic Roche model, Czechs = 1-00 kHz. 
For 0.5 Mq < M < 0.9 M^lx, the equatorial radius of Keplerian configuration of mass M, Rk(M), is, to a very 
good approximation, proportional to the radius of the non-rotating star of the same mass, Rk(M) = a R(M), with 
cins ~ «qs ~ 1.44. The value of oqs is very weakly dependent on the mass of the crust of the quark star. Both a's are 
smaller than the analytic value aR OC he = 1.5 from the relativistic Roche model. 

Key words, dense matter - equation of state - stars: neutron - stars: rotation 



1. Introduction 

Because of their strong gravity, neutron stars can be very 
rapid rotators. In view of the high stability of pulsar fre- 
quency (even the giant glitches produce relatively small 
fractional change of rotation frequency, Ss 10~ 5 ), one can 
treat pulsar rotation as rigid. The frequency / of stable 
rotation of a star of gravitational mass M and baryon 
mass lower than the maximum allowable for non-rotating 
stars is limited by the (Keplerian) frequency /k of a test 
particle co-rotating on an orbit at the stellar equator. 
The relation between /k and stellar gravitational mass 
M, /k = /k{M), depends on the (unknown) equation 
of state (EOS) at supranuclear densities. Both quanti- 
ties, M and /, are measurable, and the condition im- 
plied by a measured frequency / b s of a- pulsar of mass 
M, / bs < fn(M), could be used to constrain theoretical 
models of dense matter (for a recent review of theory of 
dense matter see lHaensel et al.l 120071 ) . Numerical calcula- 
tion of Jk(M) requires precise, time consuming 2-D cal- 
culations of stationary rotating configurations in general 
relativity. Therefore, the search for a sufficiently precise ap- 
pro ximate but universal f o rmula for f^(M) is of great inter- 
est. lLattimer &; Prakashl (|2004t ) proposed an approximate 
empirical formula f K {M) « C(M/M Q ) 1 / 2 (i?/10 km)" 3 / 2 , 
where R = R(M) is the circumferential radius of static 
star of mass M, and C = 1.04 kHz does not depend on 



the EOS. In the present paper we calculate precise 2-D 
models of rapidly rotating neutron stars and quark stars 
with different EOSs. We use the relativistic Roche model 
(|Shapiro et all I1983J) to motivate the empirical formula for 
/k(M) proposed by Lattimer fc Prakashl (|2004D . We calcu- 
late the optimal value of prefactor C and we establish limits 
to the validity of the empirical formula. 

As of this writing, the maximum ro tation frequency of a 
pulsar is 716 Hz (P SR J1748-2446ad. lHessels et~alll2006h . 
iKaaret et alj (|2007f ) reported a discovery of oscillation fre- 
quency 1122 Hz in an X-ray burst from the X-ray transient, 
XTE J1739-285, and concluded that "this oscillation fre- 
quency suggests that XTE J1739-285 contains the fastest 
rotating neutron star yet found" , but this observation has 
not been confirmed or reproduced. 

The problem of the constraint f^r( M) > / bs was al- 
ready considered bv lShapiro et alj (jl983l ) after the epochal 
discovery of the first millisecond pulsar P SR 1937+214 
with f b s = 641Hz (jBacker et all |1982[ ). IShapiro et all 
(|l983l ) used a formula for f max (M) based on the relativistic 
Roche model. After the announcement of ill-fated discov- 
ery of a 2 kHz pulsar, this formula was used to show that 
nearly all EOSs of dense matt er existing at that ti me were 
ruled out by this observation (jShapiro et al.l . fl989l ). 

In Sect. [J] we summarize results obtained with the rel- 
ativistic Roche model. Sect. [3] contains description of real- 
istic EOSs of nuclear matter, whereas Sect. 2] provides the 
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assumptions and methods used to calculate rotating stellar 
models. In Scct.[5]we check the validity of empirical formula 
against results of precise 2-D calculations for ten realistic 
EOS for neutron stars. Hypothetical self-bound quark stars 
with normal crust are considered in Sect. [SI Sect. [7] presents 
the static and rotating configurations in the mass-radius 
plane. In Sect. [8] we derive approximate relations between 
the circumferential radius of a static configuration and that 
of a Keplerian configuration of the same gravitational mass, 
for neutron stars and quark stars. Discussion of our results 
is presented in Sect. M 



2. Relativistic Roche model 

There exists an instructive model of neutron stars for 
which an analytic formu la for /k(M) can be obtained 
([Shapiro et al.Hl983l . 1989). It is a relativistic Roche model, 
in which the mass of the star is assumed to be strongly 
centrally condensed. Consider a continuous sequence of sta- 
tionary configurations of constant gravitational mass M, 
and rotation frequencies ranging from zero to /k- Let the 
circumferential radius of non-rotating configuration be R. 
Under the assumption of an extre me central mass conden- 
sation, IShapiro et all (11983L Il989h found an equation sat- 
isfied by the coord inates of the stellar surface (Eq. (2) of 
IShapiro et al.l[l989l ). In the special case of the stellar equa- 
tor, this equation implies that for normal equilibrium con- 
figurations rotating uniformly at /, the equatorial circum- 
ferential radius i? e q satisfies 



2GM 



2GM 
R 



(1) 



where the extreme central concentration of matter implies 
the gravitational mass can be treated as constant and equal 
to the static value M. The left-hand-side of Eq. JT]) reaches 
a minimum at R cq = (GM/Att 2 f 2 ) 1 ^ 3 . For the solution to 
exist for a given /, the value of the left-hand-side at this 
minimum sh ould not exceed 2GM I R, w hich implies a con- 
dition on / ([Shapiro et all . 11981119891 ). 



/ < /k = 
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/ GM\ 
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Therefore, the Keplerian frequency is 



/£ och0 (M) = 1.00 kHz 



M 



1/2 



R 



10 km 



-3/2 



(2) 



(3) 



As stated in IShapiro et alj (|1989| ), "the Relativistic Roche 
model provides a surprisingly accurate estimate of the max- 
imum rotation rate along constant-rest mass sequences " for 
many realistic EOSs. 

It is easy to show that an additional relation be- 
tween Keplerian and static configuration can be obtained. 
Namely, using Eq. (Q]) , one obtains a formula expressing i? oq 
for Keplerian configuration, _Rk, in terms of R for static 
configuration of the same mass M, 



R K (M) = -R(M) . 



(4) 



The formula for f K (M), Eq. and that for R K (M), 
imply that /k(M) is equal to the orbital frequency of a test 



Table 1. Equations of state of neutron star core. N - nu- 
clcons and lcptons. NH - nuclcons, hypcrons, and leptons. 
Exotic states of hadronic matter are indicated explicitly. 
Maximum allowable mass for non-rotating stars, , and 
the circumferential radius of non-rotating stars of 1.4 M@, 
i?i.4, are given in last two columns, respectively. 



EOS 


model 


ref. 


n rstat 
max 

[M ] 


RlA 

[km] 


FPS 


N, energy 
density functional 


a 


1.800 


10.85 


GN3 


N, relativistic 
mean field 


b 


2.134 


14.22 


DH 


N, energy 
density functional 


c 


2.048 


11.69 


WFF1 


N, variational 
theory 


d 


2.136 


10.47 


APR 


N, variational 
theory 


d' 


2.212 


11.42 


BGN1H1 


NH, energy 
density functional 


e 


1.630 


12.90 


BBB 


N, Brueckner theory 


f 


1.920 


11.13 


GNH3 


NH, relativistic 
mean field 


g 


1.964 


14.20 


GMGS-Km 


N + mixed 
N-kaon condensed 


h 


1.422 


9.95 


GMGS-Kp 


N + pure 
kaon condensed 


i 


1.420 


13.20 



Re ferences for th e EOS : a -IPandharip ande &i Ra vcnhall ( 1989f ): 
b - iGlendenning] (|1985|) ; c - iDouchin fc Haensell (|200ll ); d - A14 
Argonne NN potential and Urb ana VII three body NNN poten- 
tial, from lWiringa et al.1 ([1988D ; d' - A18 Ar gonne NN potential 
with relativistic corre ctions and Urbana m odi fied UIX* NNN 
oten tial model, from lAkmal et al.l (|1998| ) ; e - iBalberg fc Gall 
19971 ); f - Paris two-body NN potentia l and U rba na UIX three 
body NNN potential, from lBaldo et al.1 l| 19971 ); g - iGlendennTngj 
( 19 81); h, i - kao n condensate models with Ux 11 
iPons et al.1 (|2000D 



-130 MeV, 



particle orbiting at r = i?K in the Schwarzschild space- 
time around a point mass M at r = 0. An approximate 
equality /k(M) ps (M, R K ) was shown to be valid 

within a few percent for normal neutron stars and quark 
stars (jBeiger et al.l . l2007t) . This relation holds strictly for 
the relativistic Roche model, 



_L ( GM \ 



1/2 



(5) 



3. Realistic EOSs of hadronic matter 



In view of a high degree of our ignorance concerning the 
EOS of dense hadronic matter at supranuclear densities 
(p > 3 x 10 14 g cm -3 ), it is common to consider a set 
of EOSs b ased on different de nse matter theories (for a 
review, see lHaensel et al.l 120071) . We used ten theoretical 
EOSs. These EOSs are listed in Tab. Q] where the basic 
informations (label of an EOS, theory of dense matter, ref- 
erence to the original paper) are also collected. 

Six EOSs are based on realistic models involving only 
nucleons (FPS, BBB, DH, WFF1, APR, GN3). Four re- 
maining EOSs are softened at high density either by the 
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appearance of hyperons (GNH3, BGN1H1), or a phase tran- 
sition to a kaon-condensed state (GMGS-Km, GMGS-Kp). 

For GMGS-Km and GMGS-Kp models, the hadronic 
Lagrangian is the same. However, to get GMGS-Kp one 
assumes that the phase transition takes place between two 
pure phases and is accompanied by a density jump, cal- 
culated using the Maxwell construction. The GMGS-Km 
EOS is obtained assuming that the transition occurs via a 
mixed state of two phases (Gibbs construction). A mixed 
state is energetically preferred when the surface tension be- 
tween the two phases is below a certain critical value. As 
the value of the surface tension is very uncertain, we con- 
sidered both cases. 

In all EOSs models (except FPS) , the core EOS was 
joine d with the DH EOS of the crust (jDouchin fc Haensell . 
1200 It ). For the FPS model, the FPS EOS of the core was 
supple mented with the FPS crust EOS of lLorenz et al.l 
(|1993l) . 

Our set of EOSs includes very different types of mod- 
els. This is reflected by a large scatter of the maximum 
allowable masses of non- rotating stars, 1.42M© ^ M max ^ 
2.21M©, and a range of circumferential radii of non-rotating 



stars with M 
!)■ 



1.4M Q , 9.95km <, R 1A < 14.22km (Table 




M [M s ] 

Fig. 2. (Color online) Same as Fig. □ but for BBB, WFF1, 
and DH EOSs. 



1.4 



1 APR 

2 GN3 

3 FPS 




M [MJ 



Fig. 1. (Color online) Precise values of Keplerian frequency 
/k (solid line) and those calculated using Eq. © (dashed 
line), assuming C = Cns = 1-08 kHz, versus stellar mass 
M. We consider masses 0.5 Mq < M < 
APR, GN3, and FPS EOSs. 



0.9 M££ and 



T 



1 BGM111 

2 GMGS-Km, GMGS-Kp 

3 GNH3 




Fig. 3. (Color online) Same as Fig. Q] but for BGN1H1, 
GNH3, GMGS-Km and GMGS-Kp EOSs. Notice that due 
to a very strong softening by the kaon condensate and si- 
multaneous constraint M < 0.9 M^j£, the GMGS-Km and 
GMGS-Kp curves do not contain kaon-condensed segments. 
Therefore, the curves for both these EOSs coincide. 



4. Calculating stationary rotating configurations 

The stationary configurations of rigidly rotating neutron 
stars have been computed in full general relativity by 
solving the Ein stein equations for stat i onary axisymmetric 
spacetime (see iBonazzola et al.l 119931 : iGourgoulhon et all 



1 19991 for the complete set of partial differential equa- 
tions to be integrated). The numerical computations 
have been performed using the rotstar code fr om the 
LORENE library ( http :/ /www.lorene.obsprmfr| ) . The 
code implement s a multi-dom ai n sp ectral method 
introduced in IBonazzola et al.l (|1998l ). A descrip- 
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1.4 



0.8 



0.6 



1 SQM1 

2 SQM2 

3 SQM3 




0.5 



1.5 



M [M s ] 



Fig. 4. (Color online) Same as Fig. \T\ but for three EOSs 
of quark stars. Quarks stars possess maximal crust, with 
bottom density pb = Pnd- To get dash lines, we used Eq. 
© with C = C QS = 1.15. 



tion o f the code can be found in iGourgoulhon et all 
(fl999h . The accuracy of the calculations has been 
checked by evaluatio n of the GRV2 and GRV 3 virial 
error indicators IGourgoulhon fc Bonazzolal 11994 
iBonazzola fe GourgoulhonI I1994D . which showed val- 
ues lower than ~ 10~ 5 . 



5. Maximum frequencies for realistic EOSs of 
neutron stars 

In Figs.[T]-l3]we compare precisely calculated Keplerian fre- 
quencies with those given by the empirical formula 



MM) « C 



M 



1/2 



R 



10 km 



-3/2 



(6) 



where M is the gravitational mass of rotating star and R is 
the radius of the non-rotating star of mass M , R = R(M). 
The optimal value of the C prcfactor is Cns = 1-08 kHz. 
The precision of the empirical formula for /k stays remark- 
ably high for 0.5 M Q < M < 0.9M^. Relative deviations 
are typically within 2%, with largest deviations of at most 
6% for the highest masses. 



6. Quark stars 

The case of strange stars, built of self-bound quark matter, 
is different from that of ordinary neutron stars. Matter dis- 
tribution within quark stars has very low density contrast 
between the quark core edge and its center. We considered 
three EOSs of self-bound quark matter, based on the MIT 
Bag Model (jFarhi fe Jaffd . ll984t IZdunikl . l2000fl . Model pa- 
rameters are given in Table [2] Quark stars are likely to have 
a thin normal matter crust, with bottom density, f>b, not 




Fig. 5. (Color online) Gravitational mass, M, versus equa- 
torial radius, R oq , for static and rigidly rotating neutron 
stars, based on the DH EOS. Solid line S: static mod- 
els (i.e., R(M)). Solid line K: Keplerian (mass-shedding) 
configurations (i.e., i?x(M)). The area, bounded by the 
S, K curves and a dash line S max K max , consists of 
points corresponding to stationary rotating configurations. 
Configurations belonging to a shaded triangular area above 
the dot line S max — K^ nax have baryon mass Mb larger than 
the maximum allowable baryon mass for non-rotating stars, 
M^. Three lines corresponding to neutron stars rotat- 
ing stably at / = 641 Hz, 716 Hz, and 1122 Hz, are labeled 
with rotation frequencies. The nearly horizontal dot line 
Si. 4 — Ki.4 corresponds to configurations with fixed baryon 
number equal to that of the non-rotating star of gravita- 
tional mass M = 1.4 M Q . Generally, solid lines connecting 
filled circle with filled square correspond to M = const., 
while dot lines connecting filled circle with open square cor- 
respond to stars with Mb = const.. For further explanations 
see the text. 



exceeding the neutron drip density /?nd ~ 4 x 10 11 g cm 3 . 
Maximum mass of normal crusts is reached for p^, = Pnd- 

First we consider quark stars with a maximum crust. 
As we see in Fig. 01 precision of empirical formula within 
the mass range 0.5 M Q < M < 0.9 M^* is as high as for 
neutron stars (typical relative deviation within 2%, largest 
deviation of about 4% at highest masses). However, the 
value of C is larger than for neutron stars, Cqs = 1.15. 

Let us consider now quark star models with less 
massive crusts. These models were constructed assuming 
P\> < Pnd- The effect on the optimum value Cqs ~ 
/ K (M) (M/M )" 1/2 (R/10 km) 3/2 turned out to be very 
small. At a fixed M, decrease of pb leads to an increase 
of /k (more compact star). Simultaneously, however, the 
static value of R(M) decreases, and therefore both effects 
cancel out to a large extent. Consequently, Cqs depends 
rather weakly on the crust mass, and in principle one may 
use Cqs = 1.15 for any crust. 
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Table 2. Parameters of the bag models for quark stars. B 
- MIT bag constant, m s - strange quark mass. For all mod- 
els the QCD coupling constant equals a s = 0.2. Maximum 
allowable mass for non-rotating stars, M^*, and the cir- 



cumferential radius of non-rotating stars of 1.4 Mq, 
are given in last two columns, respectively. 



Rl.4 



EOS 


B 

[MeV fm" 3 ] 


m a c 2 
[MeV] 


71 TStat 

™raax 

[Mq] 


Rl.4 

[km] 


SQM1 


56 


200 


1.90 


11.27 


SQM2 


45 


185 


2.02 


11.86 


SQM3 


67 


205 


1.65 


9.94 



7. Static and rotating configurations in the 
mass-radius plane 

The formulae for /k(M) are based on a one-to-one corre- 
spondence between a static configuration S, belonging to 
static boundary S of the region of rotating configurations, 
and the rotation frequency of a Keplerian configuration K 
on the K boundary. This correspondence is visualized in 
Fig. based on the numerical results obtained for the DH 
EOS. The frequency of rotation of a Keplerian configuration 
K is obtained via the mass and radius of a static configu- 
ration S with same M. Both configurations are connected 
by a horizontal line in the R cq — M plane. 

The empirical formula for the absolute upper bound 
on / o f stably rotating c o nfigur a tions for a give n EOS 
&. s (lHaensel fc Zdunikl . Il989t iFriedman et all . 11989 ; 



Sha piro et al. . 19891 : iLattimer et all 19901 ; Haensel et al 



19951) . is of a different character. It results from an (ap 



proximate but precise) one-to-one correspondence between 
the parameters of two extremal configurations, static S max 
and Keplerian K max (filled circles), and reads 



EOS 
max 



M s 



Mr. 



1/2 



pstat \ 

10 km J 



(7) 



where C is to a very good approximation independent of 
the E OS. We have C NS « Cqs = 1-22 kHz (jHaensel et all . 
Il995f ). This value is noticeably higher than Cns or Cqs, 
which determine / K (M) for 0.5M Q < M < 0.9M^. 

The functional form of Eq. Q is, in fact, exact in gen- 
eral relativity for uniform rot ation of stars with t he so- 
called minimum period EOS of lKoranda et all (|l997t ). This 
EOS contains the single parameter, e c , which is the tran- 
sition energy density between the low-density EOS with 
P = and the high-density EOS with P = e - e c . The 
value of C for the maximum mass case is 1.35 kHz. 

In Fig. [5] we also displayed the correspondence between 
the stellar configurations of the same baryon number. The 
line S max — ► K max separates " supramassive" configurations 
from the " normal" ones, which can be reached by spinning 
up a non-rotating star. The maximum rotational frequency 
for the "normal" sequences ( r eached at point K max ) has 
been discussed bv lCook et al.l (|1994al TbT) for polytropic and 
realistic EOSs ( their Tables 3 and 7 respectively). It should 
be noted that this value cannot be estimated using our 
formula for /k(M), Eq. ©, because our formula is valid 
within a restricted mass range 0.5M© < M < 0.9M^^. 

Formula (6) connects configurations of the same grav- 
itational mass M. For example, it connects M = 1.4 Mq 



non-rotating star (filled circle - Si. 4) and Keplerian con- 
figuration (filled square - K1.4), in Fig. 5. They are joined 
by a solid horizontal line. At a fixed baryon mass, Mb, M 
increases with increasing rotation frequency. Dot line con- 
necting filled circle (Si. 4) and an open square near K1.4 
contains configurations with fixed Mb, equal to that of 
a non-rotating star with 1.4 Mq. Deviation of solid line 
from the dot one visualizes the rotational increase of M at 
a fixed Mb. For a 1.4 Mq star, the fractional increase is 
equal to 1.7%, and for the maximum static mass (dot line 
Smax -* K° max ) it reaches 3.8%. 

8. Relation between R(M) and R K (M) 

Let us consider a family (sequence) of stationary configu- 
rations rotating stably at a frequency /. They form a curve 
in the M — R cq plane (see examples in Fig. [5]) . The curve is 
bound at R oq = i? m in(/) by the axisymmetric instability, 
implying star collapse into a Kerr black hole. The largest 
circumferential radius is re ached for the Keple rian configu- 
ration, i?max(/) = ifo(M). iBeieer et al.l(|2007t) have shown 
that i?max(/) is (within 2%) equal to the radius of an orbit 
of a point particle moving in the Schwarzschild space-time 
around a point (or a spherical) mass M. This implies 



-Rmax(/) 



(- 



GM 



1/3 



For convenience we introduce a frequency /o, 



= J_ / GMq 
10 27rV(10km) 3 



= 1.8335 kHz 



(8) 



(9) 



Validity of the empirical formula, Eq. ([S]), suggests then an 
approximate proportionality 



Rk(M) w a R(M) , a = ( £ 

Jo 



-2/3 



(10) 



For the extreme relativistic Roche model, Rk(M) is strictly 
proportional to R(M), and aR, oc h c = 1.5, Eq. For neu- 
tron stars and quark stars with crusts, with masses within 
0.5M Q < M < 0.9M^ X , the proportionality holds within 
a few percent, as shown on Fig. [TJJ However, the best-fit 
proportionality factors are smaller than 1.5 of the Roche 
model, a NS w a Q s 1.44. 

The dependence of <iqs on the crust mass, M cr , is very 
weak. This can be explained via the effects of M cr on R(M) 
and /k(M). These effects oppose themselves: at fixed M, 
R(M) increases, and /k(M) decreases, with an increasing 
M cr . The cancellation of both effects results in an only slight 
de crease of aps w ith in crease of M cr (see Fig. [6]) . 

lLasota et al.l (|1996| ) derived an approximate relation 
between equatorial radius of a maximally rotating con- 
figuration, RJ^ and the radius of non-rotating neutron 
star with maximum allowable mass, i?^'/* . Maximally 
rotating configuration, stable both with respect to mass 
shedding and axisymmetric perturbations, is actually very 
close to that with largest mass, M™* x (in Fig. [5] they 
are indis tinguishable) . The approximate proportionality 
found by lLasota et al.1 (|l996f ) for neutron stars is -R/°\ x ~ 
1.32 i?^ f at . This relation connects two extremal config- 



urations. They have different masses, related by M, 



rot 
max 



6 
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Fig. 6. (Color online) Equatorial circumferential radius of 
the Keplerian configuration, Rk{M), vs. circumferential ra- 
dius of the static configuration R of the same gravitational 
mass M, for 0.5M© < M < 0.9M^*, for neutron stars 
(solid lines) and quark stars with crust (dash lines). Three 
straight cyan lines correspond: upper line to aR OC h c = 1-5, 
middle line to a = ons = 1-44, and bottom line a = 1.39. 
Neutron stars: Color of a curve for a given EOS coincides 
with that of the EOS label (APR,. . . ,WFF1). Quark stars: 
nearly straight green, blue, and red solid lines, located in 
the lower bundle, correspond to the SQM1, SQM2, and 
SQM3 EOSs of Table[2]with a maximum solid crust. Green 
dash lines in this bundle describe results obtained for the 
SQM1 EOS of quark core and low-mass crusts: 0.3 M crjmax 
(middle green line) and 0.06 M cr max (upper green line). 



1.1 8M^* (jLasota et al.l . ll996D . In contrast, Eq. (TO]) con- 
nects normal configurations of neutron stars and quarks 
stars with same gravitational mass and holds for 0.5M Q < 
M < 0.9Af± a l 



9. Discussion and conclusions 

We have tested empirical formula for Keplerian (mass shed- 
ding) frequency of neutron star of mass M, proposed by 
lLattimer fc Prakashl (|2004[ ). Using numerical results of pre- 
cise 2-D calculations, performed for ten representative re- 
alistic EOSs of dense matter based on different dense mat- 
ter models, we find prefactor Cns = 1.08 kHz, slightly 
higher than 1.04 kHz proposed by lLattimer fc Prakashl 
(120041) . With our prefactor, the formula is quite precise 
for O.5M < M < 0.9M max (typically within 2%, maxi- 
mum deviation occurring for highest M not exceeding 6%). 
Quark stars can reach larger /k(M) than neutron stars. 
With a maximum crust on quark stars, we get Cqs = 
1.15 kHz. The value of Cqs docs not depend significantly on 
the crust mass, and can be used also for bare quark stars. 
We notice that both Cns an d Cqs are significantly larger 



than for the relativis tic Roche model, CR OC he = 1-00 kHz 
(jShapiro et al.l . ll983D . 

Using an approximate but quite p recise Schwar z schild - 
like formula, relating M, i?K, and /k (|Beiger et al.l . l2007t ). 
we show that to a very good approximation the mass- 
shedding radius at a given M is proportional to the static 
radius R{M), provided 0.5 M Q < M < 0.9 M^*. For neu- 
tron stars and quark stars we obtain the best-fit propor- 
tionality factor ons ~ «qs ~ 1-44. These proportionality 
factors are smaller than the exact factor 1.5 obtained for 
the relativistic Roche model. 

Concluding, we derived a set of empirical formulae, 
expressing Keplerian frequency and equatorial radius of 
Keplerian configuration in terms of the mass and radius 
of normal configuration of the same mass. These formu- 
lae can be used for masses 0.5M Q < M < 0.9M max . The 
formulae are approximate but quite precise, and therefore 
might be useful for constraining the EOS of dense matter 
by the observations of pulsars. 
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